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Abstract 

The non-relativistic approximation of the quark-meson-coupling model has 
been discussed and compared with the Skyrme-Hartree-Fock model which in- 
cludes spin exchanges. Calculations show that the spin-exchange interaction has 
important effect on the descriptions of finite nuclei and nuclear matter through 
the Fock exchange. Also in the quark-meson-coupling model, it is the Fock 
exchange that leads to a nonlinear density-dependent isovector channel and 
changes the density-dependent behavior of the symmetry energy. 
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1. Introduction 

Both relativistic and non-relativistic energy density functional models have 
been developed to describe the many-body system of atomic nucleus [l| . Accord- 
ing to the Kohn-Sham theorem [2j], models leading to the same energy functional 
are equivalent as far as physical results are concerned [3j. Indeed, systematical 
calculations have shown that both relativistic and non-relativistic models give 
quite similar results for the ground-state properties and low-energy dynamics 
of stable nuclei However, significant differences are found when calculations 
go to the extreme cases of high isospin and high density [l|, [B| . The origin of 
these differences is still not clear. 

In the non-relativistic approach, the Skyrme-Hartree-Fock (SHF) model has 
been extensively employed in the investigation of nuclear structure [H, 0| since 
the first work of Vautherin and Brink [6]. The Skyrme potential contains ten 
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parameters, but only six of the parameters have been well determined Q- This 
leads to many sets of parameters with different emphases on the descriptions of 
nuclear properties d, Q . It was tested by Stone et al. @ that about one third 
of existing Skyrme parameter sets give the monotonic increase of the symmetry 
energy with density while other two third give the decrease of the symmetry 
energy after certain supra-saturation density. It should be mentioned that pre- 
dictions for the density dependence of the symmetry energy depend strongly 
on models and parameters employed 

EES EH- As examples, the Variational- 
Many-Body [12Hl4j and Dirac-Brueckner-Hartree-Fock 15 1 theories predict that 
the symmetry energy turns to decrease after certain supra-saturation density, 
while RMF 16, 17 1 and Brueckner-Hartree-Fock models present that the 
symmetry energy increases monotonically with increasing density. Therefore, it 
is encouraging to investigate possible physics origins for the different behaviors 
of the symmetry energy at high densities [l9l - 21 1 which may shed light on our 
understanding of nuclear force. 

Based on the mean-held approximation but including the response of the 
quark structure of the nucleon in the nuclear environment, the relativistic quark- 



meson-coupling model (QMC) was proposed by Guichon et al. 22, 23]. In 
Ref. [2J] an interesting similarity was found by comparing the energy density 
functionals of the SHF model and the non-relativistic approximation of QMC 
models, providing an understanding of the physical origin of the Skyrme force 
within the relativistic approach. With the restriction of N = Z, the authors 24 \ 
obtained a set of parameters of the Skyrme force that gives reasonable descrip- 
tions of nuclear matter and finite nuclei. The work suggests that investigations 
regarding the physics of effective nucleon-nucleon interactions remain interest- 
ing. In Ref. [24(, the Skyrme spin-exchange terms that involve the x±, X2, and X3 
parameters are neglected. The present discussion will include the spin-exchange 
terms. It will be shown that these terms have important effects on nuclear prop- 
erties at high isospin and high density, particularly on the density dependence 
of the symmetry energy. 



Energy density functionals of QMC and Skyrme Hartree-Fock mod- 
els 

JMC model gives a relativistic quark-level description of nuclear mat- 
I23j |. In the model, a baryon in nuclear medium is assumed to be 
a non-overlapping static spherical bag in which quarks are coupled to meson 
fields in the mean-field approximation 22h25|. For modeling, one can choose 
an isoscalar-scalar field (ex) for the medium-range attraction, an isoscalar- vector 
field (w) for the short-range repulsion, and an isovector-vetor field (p) for isospin 
channel [24|. Detailed discussions on the Hartree-Fock (HF) approximation of 
the QMC model can be found in Refs. (^llHSl- 

As has been given in Ref. for finite nuclei, the non-relativistic HF ap- 
proximation of the QMC density functional can be written as follows 
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where M is for the mass of the nucleon, t = t„ + r p for kinetic energy density, 
and p = p n + p p for density with n and p denoting the neutron and proton, 
respectively. G a = gl/m 2 a , G u = gl/m 2 ^ and G p = g 2 p /m 2 are defined with g^, 
j u , g p for coupling constants and m CT , m w , m p for the masses of the mesons. 
The parameter d is the scalar polarizability of the nucleon [24] . In the numerical 
calculations of the present work, we take d = 0.150 fm from Ref. j24j. /i s and 
ply are isoscalar and isovector magnetic moments of the nucleon, respectively. 
Ps = P P + Pn and p v = p P - Pn (p P = 2.79pN, Pn = -1.91/xn). Jn and J p 
are the spin density of the neutron and proton, respectively, with the total spin 
density of J = J n + J p . 

In Refs. @, H^] it was shown that the six elementary parameters (to, t\, 
t 3 , xq and Wo) of the Skyrme force can be understood in the framework of the 
QMC model. In the present work, we study in a more detailed energy density 
functional of the SHF model with spin-exchange parameters (x\, x%, x 3 ) which 
can produce effects through the Fock exchange. We can have a general form of 
the Skyrme energy density functional as [30 ] , 
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and 



H eS = ^[(2 + x 1 )t 1 + (2 + x 2 )t 2 ]pT 
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n so = -^W (pV • J + p„V • J„ + p p V • J p ) . (9) 



The QMC and SHF models have the same density-dependent forms in the 
effective- mass term H e g and the finite-range term Hfi n - By comparing H e s and 
"Hfin, we obtain the following relations, 
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Comparing the energy density functionals of the QMC and SHF models, 
we find that main differences in density dependence appear in the terms which 
involve the parameters t% and X3 in the zero-range term "Ho and the density- 
dependent term H3. For a qualitative study, we expand the Ho+%3 term around 
the saturation density in respect that the density range of finite nuclei is around 
the saturation density. For the coefficient of the p 2 term in the corresponding 
%o + %3 term, the first and second orders of the expansions give 
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where po is the saturation density. Similarly, the expansions of the coefficients 
of the (p n — p p ) 2 term give the following relations, 
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From Eqs. (TTSl) and (JT7J) , the parameter x 3 can be obtained immediately. 
The parameters xq, to and t 3 are related to the parameter a. We can use the 
symmetry energy of the nuclear matter at the saturation density to determine 
the parameter a. Fig. [T] plots the variation of the symmetry energy with chang- 
ing the parameter a. At each given a value, the parameters xq, to and t 3 are 
obtained using Eqs. (|T4l) -(fT6 |) . and the SHF symmetry energy |30j is calculated 
with these parameters including those given by Eqs. (TITj]) -(ri"3 |) , The QMC pa- 



rameters are taken from Ref . 24] . It is well accepted that the symmetry energy 
of the nuclear matter at the saturation density is about 30 MeV, which leads 
to an a value of about 1/6, that is consistent with the usual value in the SHF 
model, see Fig.[T] With the a value and the QMC parameters we can determine 
the xq, t and t 3 using Eqs. (fT^j) - (fTTjj) . 

For completeness, we take the Wo parameter corresponding to the spin-orbit 
term U so as Ref. (H 
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With the QMC parameters at m a = 700 MeV [24|, all the Skyrme param- 
eters including x%, X2 and x 3 have been derived, listed in Table [1] Similar to 



the result given in Ref. [24(, the obtained values of the xq, t , ti, t 2 , t 3 and 
Wo parameters are close to the SkM* force. In Ref. [24j, the parameters x\, x 2 
and x 3 are set to be zero. However, it was pointed out that the parameter x 3 



has significant effect on the symmetry energy 31) . There already exist many 
sets of Skyrme parameters. It is not our aim to derive a new set of parameters. 
In the present work, we extend the understanding of the Skyrme force in the 
framework of the QMC model, which has been done primarily in Ref. [24[ . 



3. Calculations and discussions 

As the first step, we compare the Hartree-Fock calculations of the QMC 
and Skyrme models with different parameters for the nuclear matter, shown in 
fig. [5J It is found that the calculations give similar properties of the nuclear 
matter around the saturation density. However, different behaviors are seen in 
densities away from the saturation. The Sqmc parameters which contain the 
spin-exchange parameters xi, x 2 and £3 give a similar behavior of the symmetry 
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Table 1. The Skyrme parameters obtained in the present work (labeled by Sqmc) from the 
QMC model [sj. 



Skyrme Parameters 

t (McV fm 3 ) -2648.19 

ti (MeV fm 5 ) 371.07 

t 2 (MeV fm 5 ) -121.644 

t 3 (McV fm 4 ) 15553.495 

x 0.60146 

xx 0.2697 

x 2 -0.23701 

x 3 0.6968 

Wo (MeV fm 5 ) 104.498 

a 1/6 
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Figure 1. The symmetry energy for nuclear matter at the saturation density as a function of 
the a parameter, calculated using the SHF model with the parameters obtained through the 
QMC model. The dashed line indicates the commonly-adopted symmetry energy of a B = 30 
MeV. 
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P (fm 3 ) p (fm~ 3 ) 

Figure 2. (Color online) Comparisons between the QMC and Skyrme calculations of the 
properties of the nuclear matter, i.e., incompressibility K^, pressure P, binding energy Eg 
and symmetry energy E Bynl within the Hartree-Fock approximation. For the SHF calculations, 
SkM* , Sqmc and the Skyrme parameters given in Ref. |24| have been used. 



energy against density as the QMC model. The symmetry energy curve calcu- 
lated by the Skyrme force neglecting the x±, 22 and £3 parameters [3] is much 
more stiff than that obtained by the QMC model. A soft symmetry energy 
against density seems more acceptable from the recent analysis of experimental 



data [11 1 



As a further test, we calculate the properties of finite nuclei by the Sqmc pa- 
rameters. We use the SHF-BCS code of ev8 [32( with a self-consistent two-body 
center-of-mass correction. Calculations for the double magic nuclei have shown 
a significant improvement, compared with those in Ref. 2J] where the param- 
eters xi, X2 and X3 were set zero. We have also calculated binding energies for 
the long chain of Sn isotopes from 100 Sn to 132 Sn, for which the paring strength 
is fitted by reproducing the neutron separation energies of Sn isotopes, shown 
in Fig. [3] We see that reasonable agreements between data and calculations 
are obtained with the Sqmc parameters. Here we should remind that the Sqmc 
parameters are obtained by comparing the energy density functionals of the 
QMC and SHF models for the nuclear matter without fitting the properties of 
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Figure 3. (Color online) Difference between experimental and calculated binding energies 
for even-even Sn isotopes. The calculations are done using SHF+BCS with different sets of 
Skyrmc parameters. 
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pi gure 4. (Color online) Neutron and proton density distributions for ^^Pb (indicated by p n 
and p v , respectively) calculated with Sqmc, SLy4, SkM* and the Skyrme parameters given in 
Ref. [24| . compared with the experimental charge density |33j |. 
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finite nuclei, while other sets of the Skyrme parameters were determined by the 
fitting procedure to the experimental data of finite nuclei (including the data of 
binding energies). 

Fig. 2] shows the calculated neutron and proton density distributions for 



208 Pb, compared with the experimental charge density [33[. The Sqmc force 
gives a similar density distribution to those obtained by the SLy4 and SkM* 
forces, agreeing well with the data. We see that the calculations with the inclu- 
sion of the spin-exchange terms are improved noticeably. 

In mean-field models, the concept of the effective mass of the nucleon mainly 
describes effects related to the nonlocality of the underlying nuclear interactions 



and the Pauli exchange effects in nuclear system |34| . The effective mass for the 



nuclear matter in the SHF model can be written as [3 
K 



( h 2 l 
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where the index "/" indicates the proton or neutron. In the QMC model, the 
non-relativistic effective mass can be extracted when one collects the pr terms, 
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The first two terms are the Hartree contribution while the remaining terms come 
from the Fock exchange. 

In the SHF model, the effect from spin-exchange terms with the parameters 
xi, X2 and X3 appear explicitly in the effective mass of Eq. (fT9|) . For the QMC 
model, if the effect from the Fock exchange is neglected, the non-relativistic 
effective mass contains only the isovector term. (The isovector effective mass 
can be obtained by setting pj = in Eq. (l20|).) If the Fock term is taken 
into account, the p meson has contribution to both the isovector and isoscalar 
effective mass. The non-relativistic effective mass calculated by the QMC model 
for the symmetric nuclear matter is shown in Fig. [5] 

The longstanding problem remains about the neutron-proton effective mass 
splitting, which is crucial for extracting the density dependence of symmetry 
energy from nuclear reaction using transport models [a, l35j . The dominant 
contribution to the neutron-proton effective mass splitting is from the isovector 
channel of the energy density functional. Fig. |6] gives the changes of the effec- 
tive masses with the isospin, compared with SHF calculations. The behavior of 
effective masses in the QMC calculations is similar to that of the SkM* calcu- 
lations. It seems more acceptable recently that the effective mass of neutron is 
bigger than that of proton [3g, [37| . 
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Figure 5. (Color online) The density dependence of the non-relativistic effective mass, m* /m, 
calculated by the QMC model for the symmetric nuclear matter. 
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Figure 6. (Color online) Calculated non-relativistic effective masses (m* /m) of the proton and 
neutron against the isospin in the QMC calculations, compared with the SHF calculations with 
Sqmc, SLy4, SkM* and the Skyrme parameters given in Ref. |24l . The calculations are done 
at the saturation density of p = 0.16 fm -3 . 
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In Ref. [38[ it was shown that the three-body term has an important influence 
on the symmetry energy curves. In Fig. [2 we have already given the symmetry 
energy for the nuclear matter, calculated by the QMC and SHF models within 
the Hartree-Fock approximation. In order to gain a more detailed insight into 
the nature of the symmetry energy and the effect from the Fock exchange, 
we make further analysis. The QMC symmetry energy for the nuclear matter 
within the Hartree approximation is written as, 

<™ = 6M ( — ^ + 6M^ ( — ^TTdGTp + T^' (21) 

If the Fock exchange is taken into account, the QMC Hartree-Fock symmetry 
energy for the nuclear matter becomes, 

HF h 2 /3vr 2 V a ifZir 2 ^ [G^ G a 
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In the SHF model, the symmetry energy for the nuclear matter has the following 
form [H, 
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The symmetry energy contains two contributions: 1) from the kinetic en- 
ergy; 2) from the interaction. It is seen that the QMC and SHF models have 
similar density-dependent forms for the kinetic energy part which originates 
from different Fermi momenta of nucleons and the momentum dependence of 
the isoscalar potential. The momentum dependence of the isoscalar potential 
has been well analyzed with experimental data [39l . [io| . The different behav- 
ior of the density dependence in the symmetry energy would be caused by the 
interaction part which is from the isovector potential of the energy density func- 
tional. From Eq. (|21[) in the Hartree approximation, we see that the symmetry 



energy contributed from the interaction part, -£p, has a linear form respecting 
the density which dominates the symmetry energy. The Fock exchange leads 
to a nonlinear density-dependent interaction part for the symmetry energy, see 
Eq. ([22|) for the QMC model. Fig. [7] displays the curves of the symmetry en- 
ergy with and without the Fock exchange. It is seen that the inclusion of the 
Fock exchange gives a strong density-dependent symmetry energy. In the fig- 
ure, we have also given the symmetry energy curve of the Hartree calculation 
with re-fitted parameters (readjust the QMC coupling constants by fitting the 
properties of the nuclear matter at the saturation) to see if the effect of the Fock 
exchange can be included by adjusting parameters. We see that the effect from 
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Figure 7. (Color online) Symmetry energy curves calculated by the QMC model in the Hartree 
and Hartree-Fock approximations. The dashed line is obtained by the Hartree approximation 
using the same parameter values as in the Hartree-Fock calculation. The dotted line is also for 
the Hartree calculation but with parameters readjusted by re- fitting nuclear matter properties 
at the saturation density. 



14 




15 



the Fock exchange can not be included by readjusting parameters. It has been 
pointed out that models within the Hartree approximation give nearly-linear 
symmetry energy [l6j |. 

It is seen from Eq. (|2"2"|) that the a meson brings a density dependent contribu- 
tion to the interaction part of the symmetry energy through the Fock exchange. 
In Fig. [5J we show the contributions of the mesons to the symmetry energy 
in the QMC-Hartree and QMC-Hartree-Fock models. We see that, with the 
Fock exchange included, the contributions of the p and cu mesons have nearly 
linear but opposite trends against the density. These two parts can be canceled 
partly. The a meson through the Fock exchange produces a strong density de- 
pendence, which increases first with density and then decreases drastically in 
higher density, leading to a soft behavior of the symmetry energy. 

We further study effects on the symmetry energy from the density-dependent 
meson-nucleon exchange interaction. If the QMC model is compared with RMF 
model, one would find that the main difference is from the scalar field 4l| . The 
scalar field in the QMC model has a nonlinear density dependence, which is 
derived by considering subnucleonic freedoms 2J, |25| . 141] . This nonlinear be- 
havior originates from the scalar polarizability parameter d of the nucleon. In 
the MIT bag model, it can be well represented by the bag radius 24, 28, 29j |. 
Fig- M plots the density dependence of the symmetry energy with assuming the 
different values of the parameter d. The standard value of the parameter d is 
0.150 fm which corresponds to a bag radius of 0.8 fm 24j. If we set d — which 
is equivalent to assuming a point particle for the nucleons, the meson exchange 
interaction becomes density independent, which gives approximately linear sym- 
metry energy, see Fig. [9] With increasing the parameter d, the symmetry energy 
becomes soft against the density. However, the Hartree calculation still give stiff 
linear symmetry energy. By comparing these calculations, we can conclude that 
the nonlinear scalar field can soften the symmetry energy curve through the 
Fock exchange interaction. As mentioned already in the introduction of the 
paper, the symmetry energy is still a debating problem. Different models give 
quite different predictions of the symmetry energy However, a 

soft symmetry energy seems to have been supported by the recent analysis of 



experimental data [11 



4. Summary and conclusion 

With the inclusion of spin-exchange in the Skyrme force, we have discussed 
the relation between the Skyrme and meson-nucleon exchange forces. As in 
Ref. [24]] . the Skyrme force can be understood within the QMC model. By 
comparing the energy density functionals, a set of Skyrme parameters including 
the Skyrme spin-exchange parameters Xi, x-i and x% is obtained, which can 
give reasonable results for the properties of nuclear matter and finite nuclei. 
However, the determination of the parameters is not the aim of the paper. 
Many sets of the parameters for the Skyrme force have existed with fitting to 
the experimental data of finite nuclei also. 
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Figure 9. (Color online) Symmetry energy curves calculated by the QMC model in the Hartree 
and Hartree-Fock approximations with assuming different values of the parameter d. The 
dashed lines are obtained by the Hartree approximation and the solid lines are obtained by 
the Hartree-Fock calculation. All the calculations use the same parameters as QMC Hartree- 
Fock except the parameter d. 
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We have paid special attention to the isovector channel of the nuclear energy 
density functionals. The QMC and SHF models have a similar density depen- 
dence for the isoscalar channel, but have different forms for the isovector chan- 
nel. It is found that the Fock exchange leads to a nonlinear density-dependent 
isovector channel in the QMC model, which makes a soft symmetry energy. The 
nonlinear density dependence is produced by the scalar field which is derived 
by considering subnucleonic freedoms. 
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